SCATTERING PHASE ASYMPTOTICS 
WITH FRACTAL REMAINDERS 



SEMYON DYATLOV AND COLIN GUILLARMOU 

Abstract. For a Riemannian manifold (M,g) which is isometric to the Euclidean 
space outside of a compact set, and whose trapped set has Liouville measure zero, 
we prove Weyl type asymptotics for the scattering phase with remainder depending 
on the classical escape rate and the maximal expansion rate. For Axiom A geodesic 
flows, this gives a polynomial improvement over the known remainders. We also 
show that the remainder can be bounded above by the number of resonances in 
some neighbourhoods of the real axis, and provide similar asymptotics for hyperbolic 
quotients using the Selberg zeta function. 

In this paper, we derive high energy asymptotics for the spectral shift function, or 
the scattering phase, with remainders which depend on the dynamic behavior of the 
underlying classical system. Our main setting is an n + 1 dimensional Riemannian 
manifold (M, g) which is Euclidean near infinity, but we shall also consider in the 
second part of the paper the case of convex co-compact hyperbolic quotients T\EI n+1 
using a different approach. If A is the (nonnegative) Laplace-Beltrami operator on M 
and A is a reference operator (such as the Laplacian on the free space when (M, g) is a 
metric perturbation of the Euclidean space), then the spectral shift function is the trace 
of the difference between the spectral projectors of A and A - see (1.1) and (1.2). The 
scattering phase is the phase of the determinant of the relative scattering matrix. In a 
very general setting, these two functions are equal almost everywhere on the absolutely 
continuous spectrum, as shown in the classical work of Birman-Krein [BiKr], thus we 
will denote both of them by the symbol s(z), where z G M is the frequency. 

Our recent work [DyGu] contains a local Weyl law with a remainder expressed in 
terms of classical escape rate - in particular, if the metric has constant curvature — 1 
near the trapped set K and the Hausdorff dimension of K is dim^(fC) = 25+1, 
then the remainder is 0(z s+ ). If M is Euclidean near infinity, we use the commutator 
method of Robert [Ro] to deduce an asymptotic expansion of s(z) with the same 
remainder - see Theorem 1. This seems to be the first result on spectral asymptotics 
with a fractal remainder which is intermediate between the usual Weyl law and the 
non-trapping case where there is a complete expansion at high frequency. By contrast, 
the best known remainder for Weyl law on compact manifolds with chaotic geodesic 
flows is 0(z n / log z), see Berard [Be]. 
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We next give several other estimates on the spectral shift function in different set- 
tings, with remainders related to the one discussed in the previous paragraph. First of 
all, the function s(z) admits a meromorphic extension into the whole complex plane, 
and the asymptotic behavior of s(z) for large real z is related to the distribution of its 
poles, known as scattering poles or resonances. In the Euclidean near infinity setting, 
we get an asymptotic expansion of the scattering phase with remainder depending on 
the number of resonances in the ball of size clog(z) centered at z, for some c depending 
on the injectivity radius - see Theorem 2. Next, for the case of hyperbolic surfaces 
T\EI 2 , we get an asymptotic expansion of s(z) using the Selberg zeta function - see 
Theorem 3. 

1.1. Euclidean near infinity case — an application of [DyGu]. Let g and g\ be 

two smooth metrics on a manifold M. When there exists a compact set iV C M such 
that each (M \ N,gj) is isometric to R n+1 \ B(0,R ) for some _R > 0, we define (up 
to a constant) the spectral shift function s(z) for the pair (A 91 , A 50 ) by duality: for 
all <p G Cq°(R), 



where it is straightforward to see that the trace is well defined, using that g = g\ 
outside a compact set. In this setting, s(z) will be a smooth function on (0, oo) and it 
is equal, up to a multiplicative constant, to the phase of the relative scattering operator 
by the work of Birman-Krein [BiKr]. 

We shall also be interested in the case where (Mo, go) and (Ml, gi) are two Riemann- 
ian manifolds such that for each j = 0,1 there is a compact set Nj C Mj with Mj \ Nj 
isometric to IR n+1 \ B(0,R ). The spectral shift function for the pair (A 9l , A 90 ) can 
then be defined using the black-box trace [Sj97] by the formula 



where the trace Trt,b is defined in (2.1). 

Before we state our result, we need to define a few geometric quantities. Let (M,g) 
be a Riemannian manifold such that, for some compact set N C M, the end (M\N,g) 
is isometric to IR n+1 \ B(0,R ) for some R > 0. Let g l be the geodesic flow acting 
on the unit cotangent bundle S*M. The trapped set K C S*M is defined as follows: 
(m, v) lies in K if and only if the corresponding geodesic g t (m, v) lies entirely in some 
compact subset of S*M. Denote by \il the Liouville measure on S*M generated by 
the function ^Jp(m, v) = \u\ g . The set T(t) of geodesies trapped for time t > is 




(1.1) 




(1.2) 



defined as 



T(t) := {(m, v) e S*M \meN, 7i(g\m, v)) e N}. 



(1.3) 
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where ir : S*M — > M is the canonical projection. We define the maximal expansion 
rate of the geodesic flow as follows: 

A max := limsup T^log sup ||d^*(m,z/)||. (1.4) 

\t\^+oo \t\ (m,i/)ef(t) 

The norm on the right hand side is with respect to the Sasaki metric. Using our recent 
work [DyGu], we obtain 

Theorem 1. Let (M, g) be a smooth Riemannian manifold for which there is a compact 
subset N such that M\N is isometric to IR n+1 \ B(0, Rq). Assume that the trapped set 
K has Liouville measure and let s(z) be the spectral shift function associated to the 
pair (A g , Agn+i) defined by (1.2). Let A max be defined in (1.4) and Aq > A max; then 
there exist some coefficients Cj such that as z — >■ oo we have for all L G N 

L 

s(z) = J2c J z n+1 - J + O(z n f i L (T(A 1 log\z\))) +0(\z\ n ~ L ). (1.5) 

3=0 

The coefficients Cj are integrals of local Riemannian invariants, they appear in the 
small time asymptotics for the local trace of the heat kernel or Schrodinger propagator 
(e.g. see Robert [Ro]); for instance, Co = C n (Vol euc i(S(0, Rq)) — Vol g (N)) where C n is 
a universal constant depending only on n. When K = 0, we recover a full expansion 
for s(z), when K ^ 0, the remainder is not better than 0(1). 

In the proof, we need to assume that M is Euclidean at infinity for (1.5) since we 
use the commutator method of Robert [Ro] to express the scattering phase as a trace 
of the spectral projector in a compact region, and then we apply our estimate [DyGu, 
Theorem 3] on the expansion as h — > of the trace Tr(A t[ 0t x\(h 2 A g )) when A > is 
in a compact interval and A G \& (M) is any compactly supported semiclassical pseu- 
dodifferential operator. It is likely that the result extends to perturbations of M. n+1 
which are not compactly supported. 



We can now specialize to particular cases where the flow is uniformly partially hy- 
perbolic in the following sense: there exists A > and an invariant splitting of TS*M 
over K into continuous subbundles 

T Z S*M = Ef © E%, VzeK 

such that the dimensions of E u and E cs are constant on K and for all e > 0, there is 
to G R such that 

Vf G E%, \\dglv\\ > e xt \\v\\, 
Vf G Ef, \\dg\vW < e £t \\v\\. 

We can define J u , the unstable Jacobian of the flow, by 



Wz eK, Wt> t , 



J u (z) :=-d t (detdg\z)\ E] 



\t=o 
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where dg 1 : E u z — > E"^tr z \ and the determinant is defined using the Sasaki metric (to 
choose orthonormal bases in E u ). The topological pressure of a continuous function 
(p : K — y R with respect to the flow can be defined by the variational formula 

P(cp) := sup (h^g 1 ) + [ (pdfi) (1.6) 

pi(iM(K) v J ' 

where Ai(K) is the set of g*-invariant Borel probability measures and h^g 1 ) is the mea- 
sure theoretic entropy of the flow at time 1 with respect to fi. Young [Yo] proved that 
for uniformly partially hyperbolic flows, the classical escape rate \im t ^ +0 o \ log(/i£,(T(t)) 
is equal to P(J U ). Combining this with Theorem 1, we obtain 

Corollary 1.1. Let (M,g) be a Riemannian manifold satisfying the assumptions of 
Theorem 1, and Ao > A max . If the geodesic flow is uniformly partially hyperbolic on 
K ^ with negative topological pressure P(J U ) < 0, then 



n+l 

s(z) = Y,c 3 z n+1 - 3 + 0(z n+p ^ A °). (1.7) 

j=0 



In particular, if g has curvature —1 near the trapped set K and the latter has Hausdorff 
dimension 25 + 1, then 

[n+l-5] 

s(z)= c 3 z n ^ + 0{z s +). (1.8) 

3=0 

When the flow is Axiom A, one has P(J U ) < by Bowen-Ruelle [BoRu, Theorem 5]. 
We refer the reader to [DyGu, Appendix B] for more discussions about the classical 
escape rate and the Hausdorff dimension of the trapped set. 

1.2. Relation with resonances. We also compare these asymptotics to the counting 
function of resonances for the Laplacian A g in regions close to the continuous spectrum. 
By definition, resonances are poles of the meromorphic continuation of the resolvent 
R(z) = (A g — z 2 ) -1 from lm(z) < to z G C as an operator mapping L 2 omp (M) to 
Ll c (M). We show 

Theorem 2. Let (M, g) be a smooth Riemannian manifold which is isometric to IR" +1 \ 
.6(0, Rq) outside a compact set N, with n + l odd, and let 71 be the set of resonances 
of A g in {Im(z) > 0}. If to is the radius of injectivity of g, and if for some c > 
satisfying do > n + 1 the following estimate holds for z > large 

Hp G K- \z - p\ < c\og(z)} = 0(z a ), (1.9) 

for some a G [0,n], then the spectral shift function defined by (1.2) satisfies the asymp- 
totics as z —7- oo 

n+l 

s( z ) = ^2cjZ n+1 - j + 0(z a ). 

3=0 
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The proof uses some ideas of Melrose [Me88] by expressing the derivative of s(z) in 
terms of resonances through a Poisson formula. Dimassi [Di] showed that s(z) has a 
full expansion if there are no resonances in a logarithmic neighbourhood of the real 
axis; our result is (in a way) a generalization of [Di] . 

Sjostrand-Zworski [SjZw07] proved that if the geodesic flow is hyperbolic on K, then 

Hp e K] \z- P \ < c} = o(z m ) (1.10) 

for all m satisfying 2m + 1 > dim^lT), where dim^ denotes the Minkowski di- 
mension (in case of pure dimensions, m can be taken to be (dim.M(K) — l)/2). For 
semiclassical operators with analytic coefficients satisfying hyperbolicity assumptions, 
Sjostrand [Sj90, Theorem 4.6] gives estimates for numbers of resonances in boxes of 
real part of size z x l 2 and imaginary part of any size between 1 and o(z); it is feasible 
that in combination with the second microlocalization analysis of [SjZw07] to handle 
the real parts, this method gives in particular the estimate (1.9) with a any number 
such that 2a + 1 > dxm M (K). It is not clear whether an estimate of the type (1.10) 
for resonances in the logarithmic regions (1.9) holds for general C°° metrics. 



1.3. Hyperbolic quotients. In the hyperbolic near infinity setting, there is no nat- 
ural model operator to define a spectral shift function, however in even dimension 
one can use the construction of the second author [Gu]. In dimension 2, Guillope- 
Zworski [GuZw] defined also a relative scattering phase s(z) by using the Laplacian on 
the funnels with Dirichlet condition on a curve separating the funnels and a compact 
region of M. In constant curvature T\EI 2 , this curve can be taken to be the bound- 
ary of the convex core, made of disjoint closed geodesies, and s(z) is related to the 
argument of the Selberg zeta function Zp(l/2 + iz) of r for z > 0. Using Zr(X), the 
estimate on its growth in strips proved by Guillope-Lin-Zworski [GuLiZw], and some 
classical results of holomorphic functions, we can compare our microlocal result to see 
that, at least without more dynamical assumptions, our estimate seems almost sharp: 

Theorem 3. Let M = T\M. 2 be a convex co-compact hyperbolic surface and let s(z) be 
the scattering phase defined in [GuZw] using the Laplacian on the funnels with Dirichlet 
condition at the boundary of the convex core. Let 5 G (0, 1) be the Hausdorff dimension 
of the limit set of T, or equivalently 25 + 1 is the Hausdorff dimension of the trapped 
set K . Then: 

1) s(z) satisfies the asymptotics as z — )■ +oo 



s\z 



47T 



\o\{N)z 2 + R(z), \R(z)\ = 0{z 6 



R{t)dt =0(z 6 
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where N is the convex core of M. 

2) the Breit-Wigner formula holds: for all z G [T — a,T + a] with T large and a fixed 

^) = ^Vol(A0z 2 + i £ ^% + 0(T*) 
Air ii \z — or 

where TZ a +i denotes the set of resonances in the disc of radius a + 1 centered at T. 

3) If T is an infinite index sugroup of an arithmetic group derived from a quaternion 
algebra, and if S > 3/4, then for all e > 0, R(z) is not an 0(z 2S ~ 3 / 2 ~ e ) . 

The statement 3) is essentially restating a result proved by Jakobson-Naud [JaNa]. 
A result similar to Theorem 3 (except statement 3)) is also true for even dimensional 
convex co-compact hyperbolic Schottky manifolds, see Theorem 5 in Section 3.2. 

When 5 < 1/2 in dimension 2, there is actually an exact formula for s(z) as a 
converging infinite sum over closed geodesies, see (3.7). 

1.4. Previous results. The asymptotics for the scattering phase has a long history, 
we will only give references related to the settings we use in this paper. The first 
result giving an asymptotics with the leading term was announced in Buslaev [Bu] 
for a perturbation of M> d by a non-trapping obstacle, this was shown with non-sharp 
remainder in Majda-Ralston [MaRal] for the strictly convex case, and by Jensen- 
Kato [JeKa] for star-shaped domains; an asymptotics with remainder 0(z d ~ 3 ) was 
proved for non-trapping compact metric perturbations in Majda-Ralston [MaRa2]. 
Then, Petkov-Popov [PcPo] gave a full asymptotic expansion in powers of z for non- 
trapping obstacle perturbations of R rf and a sharp estimate with 0(z d ~ l ) remainder was 
shown by Melrose [Me88] for general obstacles Q in odd dimensions using the counting 
function of resonances in large disks of C. When the set of transversally reflected 
geodesies is of measure 0, [Me88] obtained a second term of the form Vo\(dQ)z d ~ l with 
remainder o{z d ^ 1 ) using a Duistermaat-Guillemin/Ivrii type argument. 

Robert [Ro] proved a very general result for (not necessarily compactly supported) 
second order perturbations of the Euclidean Laplacian on M. d with a remainder 0{z d ~ 1 ) 
and with a complete expansion in powers of t for non-trapping perturbations. In the 
black-box compact perturbation setting, Christiansen [Chr] gave an estimate with a 
sharp remainder. In the semiclassical setting, the work of Bruneau-Petkov [BrPe] 
extends the results of [Ro, Chr]. The Breit-Wigner formula relating resonances and 
scattering phase in the semi-classical setting appears in particular in [BrPe] and in 
Petkov-Zworski [PeZw]. In the setting of hyperbolic surfaces with funnels, Guillope- 
Zworski [GuZw] proved a Weyl asymptotics with remainder O(z) (in fact they dealt 
with n c cusps as well, in which case there is a second term — ^zlog(z) before O(z)). 
For convex co-compact hyperbolic manifolds r\IEP +1 , the second author [Gu] gave an 
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asymptotics with remainder 0(z n ) for the Krein spectral function, and Borthwick [Bo] 
extended it to compact perturbations. 

2. Asymptotics of the spectral shift function for compact 
perturbations of the euclidean laplacian 

2.1. Asymptotics using the trace estimate of [DyGu]. Let (M,g) be a complete 
Riemannian manifold satisfying the assumptions of Theorem 1. We put (Mx,gi) = 
(M,g) and (M ,g ) to be R n+1 with the Euclidean metric; let N ± = N C M x and 
N = B(0,R ) C M , so that (M 1 \N 1 , g^ is isometric to (M \N ,g ). We define s(z) 
by duality like in (1.2) using the black-box trace: with the notation [A/]? = A — Ai, 
the black-box trace is defined by 

Tr bb (vp(VA^) -<p(y/A^)) :=[Tr(xM^)Xj)]° 1 + [Tr((l ~ X M^)Xj)]\ 

+ [TrixMyf^M 1 ~ Xi))]i 

+ Tr([(l- X ,M V / A^)(l-X,)]?)- 

(2.1) 

Here Xj £ QTO^Kj) is a function which equals 1 near Nj, and (1 — Xi) = (1 — X2) £ 
C°°(Mj \ Nj). All the traces exist by easy arguments, see [Sj97], and in the last trace, 
we use the identification Mi \ Ni ~ M \ N . The definition agrees with (1.1) when 
Mi = M 2 . In this setting, s(z) is always an analytic function on [0,oo), this can be 
seen for instance by using the Birman-Krein formula relating s(z) with the argument 
of the determinant of the relative scattering operator (which is analytic and of modulus 
1 on the real line). 

Let us first recall Theorem 4 in [DyGu, Section 5.3]. We refer to [DiSj, Zw] or [DyGu, 
Section 3] for the definition of semiclassical pseudodifferential operators \I/*(M) and 
semiclassical quantizations Op h . 

Theorem 4. Let (M, g) be a complete Riemannian manifold satisfying the assumptions 
of Theorem 1 and A > A max with A max defined in (1.4). Let P{h) = h 2 A g and let A = 
OP/i( a ) E \1/ (M) be a compactly supported semi-classical pseudodifferential operator. 
Then there exist some differential operators Lj of order 2j on T*M , depending on the 
quantization procedure Op h; with L = 1, such that for all L 6 N, all X > 0, we have 
as h — t- 

Tr(Al [0 ,A](P(/i)) =(2^)-"-^^ J Ljadfiu + h- n O(/2 L (T(A x | \ogh\) + h L ) 

j=0 l^li<A 
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where the remainder is uniform for X in a compact interval of (0, oo) and n w is the 
symplectic volume form on T*M. 

Using Theorem 4 and the commutator method of Robert [Ro], we will show Theo- 
rem 1. 

Proof of Theorem 1. Define (Mi,gi) and (Mq, go) as m the beginning of this section 
and put 8 = M x \N X , which is identified with B c := M \N . Let Xj e C^(Mj) be the 
cutoff functions used in the definition (2.1) of the black box trace. Let A be a smooth 
first order differential operator supported in 8, which is equal to A = |(m.V+V.m) on 
the support of 1— Xi, where m G M. n+l is the Euclidean coordinate in 8 . A computation 
shows that [A 90 ,^4] = A go in B c . Therefore, since A gi = A go on 8, we get, using the 
notation [Aj]± — A$ — Ai, 

/■oo 

/ <p(z 2 )z 2 d z s(z) dz =[Tr(xA 9 M A 9 M + [ Tr ((! " X;)A fl ^(A 9 .) Xi )]i 
Jo 

+ Tr([(l- Xj )[A gv AMA gj )(l- Xj )] 1 ) 

for all if G Co°((0, oo)), and both traces make sense. Using the argument of Robert [Ro, 
Appendice A], it is direct to see by a limiting argument (using cutoff functions sup- 
ported on large balls) and the cyclicity or the trace that 

Tr([(l - Xj )[^ 9j ,AMA gj )(l - xM) = Tr(Q«MA 90 )) - Tr(Q^(A 9l )) 
with Qj := ((1 - Xj)[Agj,Xj] + [ A gj >Xj](l - Xj))A. Now we deduce that 

POO 

/ <p(z 2 )z 2 d z s(z) dz = Tr(S <p(A go )) - Tr(^(A 9l )) 
Jo 

with 

Sj ■= Qj + \.,\. + Xi(i - xj)\ 

compactly supported second order differential operators. Let Xj £ C^(Mj) such that 
Xj = 1 on suppxj, and in particular XjSj = Sj. We then obtain for z > 

d z s(z)dz = z- 2 (Tr(S dU z2 (A go )xo) ~ Tr(S 1 dn z2 (A gi )xi) 

where dH z (A g ) is the spectral measure of A g ., i.e. such that <p(A gj ) = J °° (p(z)dH z (A g .) 
for all if G Cq°(R). The traces here make sense since dU z 2(A g .) are operators with 
smooth kernels for z > (see e.g. [DyGu, Section 6.2]) and Sj and Xj are com- 
pactly supported. Let e > be small and let ip G Cq°((1 — e, 1 + e)) equal to 1 on 
[1 - e/2, 1 + e/2], and let (p G Cg°((-1, 1 - e/2)) be such that ^(A) + p(A) = 1 for 
< A < 1. Then for h :— 1/z small 

s(l/h) = / ^(/iV^sfa) + Tr hh (ip(h 2 A go ) - ip(h 2 A gi )) 
Jo 

=[Tr(Tj l^h'A^xM + Tr hb ( V (h 2 A go ) - V (h 2 A 9l )) 
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where Tj := h 2 Sjip(h 2 A g .) with vp(\) = ^(A)/A in C°°((l - e, 1 + e)). In all terms 
above, one can remove the functions Xj D Y cyclicity of the trace and using XjSj = Sj. 
Observe that Tj G ty~°°(M) by using the Helffer-Sjostrand functional calculus [DiSj, 
Chapters 8-9]. Moreover, since Sj is a compactly supported differential operator, Tj 
can be decomposed into the sum of a compactly supported semi-classical pseudodiffer- 
ential operator in ty~°°(M) and of an operator Wj G /i°°\l/ _00 (M) with Schwartz kernel 
K Wj (m, m!) compactly supported in m and decaying to infinite order as ml — > oo. In- 
deed, an operator in ty~°°(M) has a smooth Schwartz kernel which is an 0{{ d ^ m ,^ )°°) 
where d(m, m') is the Riemannian distance (this comes from integration by parts in 
the oscillating integral defining the kernel). We then have 



where || • ||t> is the trace norm. We can then apply Theorem 4 with the compactly 
supported operator Tj — Wj G \I/ _00 (M) instead of Tj and we get an expansion in powers 
of h for Tr(Ti l m (h 2 A gi )) up to e>(>-> L (T(A x | log/i|))), and a complete expansion 
for Tr(T l[o i i](/?. 2 A 9o )) since go has no trapped set. The second term with tp{h 2 A g .) 
has a full expansion by [Ro, Theoreme 2.1], this is a consequence of the functional 
calculus for h 2 A gj , either using the Helffer-Sjostrand method [DiSj, Chapters 8-9] or 
the Helffer-Robert [HcRo] approach. This achieves the proof. □ 

In fact, this argument would apply similarly to any situation where the free Lapla- 
cian A go on M. n+1 is replaced by a self-adjoint operator Pq satisfying the black-box set- 
ting [Sj97], with Pq = A 9l on functions supported outside the black-box and for which 
there are expansions of Tr(AII[o, z ](Po)) i n powers of z up to 0(z n fx^T (A^ 1 log z))) 
when A is a compactly supported pseudo-differential operator. 

2.2. Asymptotics using resonances near the real axis. For compact perturba- 
tions P of the Euclidean Laplacian (this can be metric, obstacle or potential) in odd 
dimension n + 1, the resolvent R(z) = (P — z 2 )" 1 of P extends meromorphically with 
poles of finite rank from {Im(^) < 0} to C, the poles are called resonances and form 
a discrete set 71 C {Im(z) > 0}. We associate to each resonance Zj its multiplicity 
raj G N as in [SjZw91]. The following sharp upper bounds on the counting function of 
resonances for metric perturbations have been proved by Sjostrand-Zworski [SjZw91] 
and Vodev [Vo] 



Proof of Theorem 2. Let x £ Co°(M) be a function which is equal to 1 near N, and 
define the tempered distribution on R 



|Tr(H/, l m (h 2 A g .))\ < \\W,\\ Tr = 0(h°°) 



tt{p e K; \p\ <R} = 0(R n+1 ). 



(2.2) 



u 
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z 



Figure 1. The sets of resonances used in the proof of Theorem 2. We 
mark the regions corresponding to the sums v c , q c , and a c . 1Z C is the 
region below the curve, while the shaded region corresponds to A r c ,i g(^)- 

By a straightforward computation, for all x C£°(R n+1 ), Tr(x cos(£\/A K n+i)) i s a 
distribution supported at t = and its Fourier transform is a polynomial of order n. 
Thus, using the definition (2.1) of the black-box trace, we easily see that there exists 
a polynomial P n of degree n such that for all if G C£°((0, oo)) 

1 r 00 

— (u, v?) = Tr bb (^(v^) - <p(VA)) + (P n , y?> = - / (d z s(z) - P n {z))^{z)dz. 
^ Jo 

It then suffices to study the asymptotics of f* u(w)dw as z — > oo. In [SjZw94, Theo- 
rem 2], Sjostrand-Zworski proved the following Poisson formula extending [Me83] (see 
also [Zw97, Theorem 2] for a direct proof without using Lax-Phillips theory) 

u(t)=^e*l*l, t^O. 
pell 

Moreover, u(t) is classically conormal at t — 0, this follows from writing the wave 
propagator as a (classical) Fourier Integral Operator (see e.g. [HoIII, Sections 17.4 
and 17.5]): if 9 e C™(-t + S,t - 5) for 5 > small enough, 6{t) e [0, 1] and 6(t) = 1 
in [-t + 25, t - 25], then 

oo 

F(u(t)6(t))(z) ~ J2 a k zn ~ k asz^ioo (2.3) 

fc=0 

for some coefficients (which are integrals of local invariants in the metric g). Here 
J 7 denotes the Fourier transform and to is the injectivity radius of M. Notice that, 
by (2.2), if L is large enough, t L ^ pg 7^e ip '*' is a well defined distribution on K and it 
equals t L u(t) (in fact L can be taken to be n, see [Zw97, Theorem 7]). 

Let c > and define the region (see Figure 1) 

U c := {p E K; Im(p) < 1 + clog |Re(p)|}. 
We decompose u(t)(l — 9(t)) as 

u(t)(l - 6{t)) = u c (t) + v c (t), with w c (t) = e * pl ''( 1 - 
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Both terms are tempered distributions on R. From (2.2), one has v c (t) G L P (M) for all 
p > 1 if ct > n + 1 and 5 small enough: indeed 

\v c (t)\< Yl |l-^(t)|e' Im(p)|f| <C7(l-g(t)) R- c ^ +n dR< C ^~ 9 ^ - 

i 1 C \ t\ TIj 1 

This inequality also implies \v c (t)/t\ < C(l + 1 2 )' 1 and therefore we get F{v c (t)/t) G 
n C°(R). We thus obtain 

j% c {w)dw = iF(^){z)-i J ^dteL°°(R). (2.4) 

In particular, the left-hand side of (2.4) is 0(z a ) as z — » oo. 

It remains to consider the Fourier transform of u c . Let ip G Cq°((— 2,2)) be equal 
to 1 on (—1, 1), then we are going to show that 

u c (z) = 2ixa c (z) - a c * 6(z) + 2ixq c (z) - q c * 9(z), 

pen c W 1 

and a c (z) a symbol, that is there exists V such that for all k G N, there is > such 
that for all z G R 

\d k z a c {z)\<C k {l + \z\) L '- k . 

The proof follows closely the argument of Melrose [Me88]. We have as tempered 
distribution 

H jM + e-™){z) = ^4 + 2Im( ~f 2 ) 

1 2 — p| 2 \z + p| 2 

and since the resonances are symmetric with respect to the imaginary axis, we have 

f(^ t ««)( 2 ) = E«'(S) 

and these make sense as tempered distribution on R by (2.2). Then 
\L fA ,„\ „ , fo-x a\f~\\ \^r;*\L{f-, jJ\p\\\ 2Im (p)\ 



(td z ) L (u c (z) -q c * (2tt5 - 6){z)) = ^ (zd z ) L ((l - if, (— ))r^4) * (2^o - 0)0 
and one can easily see that 

\p\\\ 2Im(p) \ 



0~c,L{ Z ) : = 

pen 



is a symbol (see [Me88] for details) if L is large enough. Integrate & c ,l(z) L times, 
this defines a classical symbol with an order increasing by L, that we denote a c (z) and 
such that (id z ) L (a c (z)) = (x c> l(z), and we therefore deduce that 

Uc(z) ~ g c * (2tt5 - 6)(z) -a c * (2n5 - 0)0) = P L (z) 
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where Pl(z) is a polynomial of order L. Taking inverse Fourier transform of the left 
hand side give a distribution vanishing near t = and therefore Pl{ z ) — necessarily. 
This shows (2.5). Since for z large 2na c (z) — o c *Q{z) = 0(z~°°), we also get 

u c (z) = q c * (2tt5 - 6)(z) + 0{z~°°). (2.6) 

An easy computation yields 

Va<6, fp^dz<, (2.7) 

J a \Z-p\ 2 

and if Re(p) + l<a<6ora<6< Re(p) — 1 

'Mj,< ; (*-°)Im(,) (28) 
|^-p| 2 - (6-Re(p))(a-Re(p)) V 7 

Define A/" C| iog(z) := |j{p G 7^ c ; k — Re(p)| < l°g( ;Z )} an d set Qc( 2 ) = fn Qc(w)dw. Then 
we obtain from (2.7) and (2.8) that for w G [0, 1] and z ^> 1 



\Qc(z) - Q c (z -w)\< c(K,i og (z) + Yl 



Im(p) 



|Re(p)-2|>log( Z ) 

Using the assumption N c ± og (z) = 0(z a ), we easily deduce for z ^> 1 and to G [0, 1] 

\Q c (z)-Q c (z-w)\ = 0(\z\ a ). 
Thus iterating and using the same argument for z < 0, we get for all z, t G R 

|Q c (z) - q c (z - w )| < c(i + HX 1 + |z| + M)°. 

We multiply by 6{t) and integrate to deduce 

\27iQ c (z)-Q c *6(z)\<C(l + \z\) a . 

It follows by (2.6) that 

f u c {w)dw <C(l + \z\) a . 
Jo 

Combining this with (2.4) and (2.3), this yields that as z — > oo 

n+l 



/ u(w)dw = V c^ n+1 - J + 0{ 

JO .-_ n 



j=0 

for some Cj G 1R. This achieves the proof. □ 

Note that the proof applies for obstacles since in that case the small time asymptotics 
of the wave trace is proved by Ivrii [Iv]. 
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3. Scattering phase asymptotics for hyperbolic quotients 

3.1. Convex co-compact groups. A convex co-compact group T of isometries of 
hyperbolic space H n+1 is a discrete group of hyperbolic transformations (their fixed 
points on B = {m G IR n+1 ; \m\ < 1} are 2 disjoint points on § n = dM) with a compact 
convex core. The limit set Ar of the group and the discontinuity set fir are defined by 

A r : = { 7 ( m ) i e^ 1 ; 7 £ r} n §" , Q r :=S n \A r (3.1) 

where the closure is taken in the closed unit ball B and m G H n+1 is any point. The 
group T acts on the convex hull of A r (with respect to hyperbolic geodesies) and 
the convex core is the quotient space. The group Y acts properly discontinuously on 
fir as a group of conformal transformation and for convex co-compact groups, the 
quotient T\flr is compact. The manifold M = r\H n+1 admits a natural smooth 
compactification given by M = r\(IEP +1 U fir), where the smooth structure is induced 
by that of the closed unit ball B compactifying HP +1 . The hyperbolic metric on M is 
asymptotically hyperbolic in the sense of Mazzeo-Melrose [MaMe]. Let i be a non- 
negative function which is a boundary defining function for dM = T\flr- We say that 
a; is a geodesic boundary defining function if |<ilogx| = 1 near dM, where the norm is 
taken with respect to the hyperbolic metric on M. Such functions exist, this is proved 
for instance by Graham [Gr, Lemma 2.1]. 

An important quantity is the Hausdorff dimension of Ar 

5:=dim H A r <n (3.2) 

which in turn is, by Patterson [Pa] and Sullivan [Su], the exponent of convergence of 
Poincare series 

m G H n+1 , e" sd(m ' 7m) < oo <=^ s>5. (3.3) 
7er 

where here and below, d(-, •) is the distance for the hyperbolic metric on H n+1 . Note 
that the series is locally uniformly bounded in m G HP +1 . The Hausdorff dimension of 
the trapped set K of the geodesic flow on T\EI n+1 is given by 25 + 1, see Sullivan [Su]. 

3.2. Krein's type spectral function. In [Gu, Theorem 1.3], the second author de- 
fined on even dimensional convex co-compact hyperbolic manifolds M = r\H n+1 a 
generalized Krein function £,(z), which is a regularized trace of the spectral projector 
ri[ n 2 / /4 in 2/4 +2 2](A), and replaces the spectral shift function in this setting. As we will 
see, it is related to the scattering phase. It is defined as follows 

£(z)= [ Tr R (dn u (A)), z>0 (3.4) 
J o 

where dH u (A) is the spectral measure of the Laplacian, i.e. fjf dH u (A) = l r2 2 ^ ,(A), 
and the regularized trace is defined using Hadamard finite part: let x be a geodesic 
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boundary function in M, then if T is an operator with smooth integral kernel T(m, m!) 
such that the restriction to the diagonal m \-t T(m, m) has an asymptotic expansion 
at DM in powers of x(m), we set 



Tr R (T) = FP £ ^ / T(m,m)dvo\ g (m 

Jx>£ 



and FP stands for finite part (ie. the coefficient of e° in the expansion as e — > 0). The 
renormalized trace should a priori depend on the choice of x but in fact Tr R (dU u (A)) 
does not depend on x, see [Gu]. It is natural to use this definition in this setting 
as there is no model operator to compare with, unlike for perturbations of Euclidean 
space. The function £,'(z) has meromorphic extension to z G C and £ is a phase of a 
regularized determinant of the scattering operator S(n/2 + iz), see [Gu, Theorem 1.2]. 
The following identity is proved in [Gu, Theorem 1.3] when n + 1 is even (T(X) is the 
Euler function) 



y / n , . x = exp - 2vk£{z) + x M -=7 — — t- / rr ^w -+\ dt + cm I 

Z r {% + iz) V T{n + 1) J T(it)T(—it) / 

(3.5) 

where c G N is a constant given by the multiplicity of n/2 as a resonance, x(^0 is the 
Euler characteristic, Zr(A) is the Selberg zeta function, defined for 3?(A) > 5 (here S 
is defined in (3.2)) 



00 



I e -AmZ( 7 ) 

Zr(A) = exp I - > > — 7^ — r 

7 m=l ' 

where 7 runs over the set C of primitive closed geodesies of X, £(7) is the length of 
7 and Gj(rn) := e~% ml ^ \ det(l — -P™)| 5 if P 7 is the Poincare linear map associated 
to the primitive periodic orbit 7 of the geodesic flow on the unit tangent bundle. The 
Selberg zeta function extends to C and has only finitely many zeros in {Re(A) > |}, all 
on the real line by Patterson-Perry [PaPe]. The zeros of Zr(A) in {Re(A) < n/2} are 
the poles of the meromorphic continuation of the resolvent -R(A) = (A — A(n — A)) -1 of 
the Laplacian from {Re(A) > n/2} to C (as a map C^(M) C°°(M)). Keeping the 
convention in Section 2.2 for resonances, we say that p G {Im(p) > 0} is a resonance 
if this is a pole of the meromorphic continuation of z — > (A — n 2 /A — z 2 ) -1 from the 
lower half plane to C. The function with Gamma factors can be written as 

r(f +i z )T(<Z-iz) f l 2 x / n ^2 , 2 



rUz)r(-iz) 



(i + , 2 )...((|-i) 2 + , 2 ),(i + 0(,-)) 



as z 00. Notice that this is an odd polynomial modulo 0(z °°). Since Zr(A) 
Zr(X), we deduce that for z > 0, 

£(z) = X {M)F n+1 {z) + -Arg(Z r (n/2 + iz)) + \c 

71 Z 
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where 2iirF n+ i (z) is the explicit function in factor of x(M) in (3.5) (this is a polynomial 
modulo 0(z~°°)), and Arg(Z r (n/2 + iz)) is defined by 

f z / Z' ( - + it)\ 
A: g (Mn/2 + i Z )) = l f^^f^)*. (3.6) 

The case 5 < n/2. The series defining Zr(f + iz) converges uniformly for zgR and 
therefore one obtains directly an explicit formula (in this case c = as there is no 
resonance at n/2) 

« z, -" (M) r(n + i)X r(if)r(-«) rf,H ' (37) 

The last term in the right hand side is 0(1) as z — > oo, and by Gauss-Bonnet, we 
notice that x{M) can be written as a constant times a regularized volume Vo1r(M) 
of M, defined by Vo\ R (M) = dvol 9 (m) (cf [PaPc, Appendix A]), corre- 

sponding to the usual Weyl type asymptotics for spectral shift functions. The classical 
dynamic only appears at order z° (with oscillations) in the asymptotic. 

The case 5 > n/2. We can not use the convergence of Zr(f + iz) anymore. We 
know by Guillope-Lin-Zworski [GuLiZw] that for Schottky groups T (in particular all 
convex co-compact groups in dimension n + 1 = 2), one has for any o"o, that there 
exists C ao > such that for all z > 1 and o > cr 

hg\Z r (a + iz)\ < C ao z s . (3.8) 
Lemma 3.1. There exists N > such that for all zel, Re(Z r (iV + iz)) > 1/2. 

Proof. From the definition of Zr{\) for Re (A) > 5, we have 

Re(Z r (A)) = exp ( - ^ — Re(T A (m)) J cos ( ^ — Im(T A (m)) 

\ m=l / \m=l 

where T\{m) = e _Am£( - 7 - ) /G ! 7 (m). Let A = iV + iz and we shall take iV large. We 
claim that there is iV large enough and C > such that for all m G N, z G K 

|T w+lz (m)| ^Ce^" 1 ^. 

where 70 is the geodesic with the smallest length. This follows directly from a lower 
bound |G 7 (m)| > cq for some Co > independent of m and 7, and the estimate 

tf{7 G £; £(7) <R} = 0(e 5R ), R->oo. 

To get the lower bound of |G 7 (m)|, it suffices to write this term explicitly: associated 
to the primitive geodesic 7 there is a conjugacy class of hyperbolic isometries in T, and 
7 is conjugated by a hyperbolic isometry to the transformation (x, y) 1— > e i ^(0 1 (x),y), 
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in the half space (x,y) G H n+1 = W 1 x R+ where 7 G SO n (R). Denoting by 
ai(j), . . . , a n (j) the eigenvalues of 7 , then 

n 

G 7 (m) = det (J - e-^Ofl = JJ (l - e"" 1 ^ ^(7)™) . 

Therefore |G 7 (m)| > (1 - e - m ^ 7o) )" > (1 - e ^ (7o) ) n where £(70) > is the length of 
the shortest geodesic on M. 

We conclude that for any e > 0, there is C > such that for all z G R 



00 _ 

~\T N+ Um)\ < Ce- N ^^ 

< ' m 



m 

m=l 



and taking large enough and e small enough, we deduce the Lemma. □ 

One can obtain an approximation of the argument of Zr(f +iz) in terms of its zeros, 
this follows from a general result on holomophic functions: 

Lemma 3.2. Assume that T is a Schottky group in dimension n+1 even with Hausdorff 
dimension of limit set 5. Then for all o > fixed, there is a constant C > depending 
on a such that for all z > large and all A in the disc 

D(n/2 + iz, a) := {A G C; |A - n/2 - iz\ < a}, 

we have the identity 

4(A) ^ 1 



ZAX) 4^ A - 

where 7Z a +i is the set of zeros of Zr(A) in D(n/2 + iz, a + 1). 



< CV (3.9) 



Proof. We apply Lemma a of [Ti, Section 3.9] with the holomorphic function Zr(X) 
in a disc D(N + iz, 2iV) where iV > is a large fixed parameter (chosen independent 
of z). Using the bound \Z T {\)\ < e Cz& of [GuLiZw] for A G D(N + iz,8N) where C 
depends on N, and the lower bound \Zr(N + iz)\ > 1/2 for all z > if N is large 
enough by Lemma 3.1, we obtain that there exists C depending on C, such that for 
all A G D(N + iz, 2N) 

where 1Zn is the set of zeros of Zr(\) in the disc D(N + iz, AN). To reduce to a sum 
over 7Z a+ i, it thus suffices to use that the number of zeros of -Zr(A) in D(N + iz, 4N) is 
bounded by C"z a for some C" depends only on (by [GuLiZw] again), and |A — s\ > 1 
if I A — n/2 — iz\ < a and s G TZn \ TZa+i- This achieves the proof. □ 



We also have an estimate on the argument of Zr(n/2 + iz) by using another Lemma 
in Titschmarch's book [Ti, Section 9.4]. 
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Lemma 3.3. Fix N > a > 0. Let f(z) be a holomorphic function in {Re(z) > 0}, 
with f(z) eR for z eR and there exists K > such that |Re(/(iV + it))\ > K for all 
t £ R. Suppose that \f(o~' + it')\ < M aj t for all 1 < t' < t and a 1 > a with o > a~o- 
Then if T is not the imaginary part of a zero of f , we have 

|Arg(/(cr + iT))\ < C(\og(M^ T+2 ) - \og(K)) + 3vr/2 

for some C > depending only on o~q and N. 

We apply this to the function Zr(z) by using Lemma 3.1 and (3.8), we obtain directly 

Corollary 3.4. If T is a convex co-compact Schottky group of orientation preserving 
isometries o/H n+1 with n + 1 is even, and Zr(X) the Selberg zeta function for T, then 
for z > 1 large 

Axg(Z T (n/2 + iz)) = 0(z 5 ) 
where 5 > is the Hausdorff dimension of the limit set of V. 

In fact an integrated estimate is also true, by using a lemma of Hadamard: 

Lemma 3.5. Let T be a convex co-compact Schottky group of orientation preserving 
isometries ofM n+1 with n + 1 is even, and Zt(X) the Selberg zeta function for F and 
5 > the Hausdorff dimension of the limit set. Then for T > 1 large 

r-T 



J Arg(Z r (n/2 + iz))dt = 0(T S ). 



Proof. Let T > 1. We apply the Hadamard Lemma to the holomorphic function 
Zr(n/2 — iz) in the rectangle z £ + i[0, N], with no zeros in that rectangle: we 
get 



pN pN 

/ \og\Z T (n/2 + u-i)\du- / log \Z r (n/2 + u - iT)\du = 
Jo Jo 

Arg(Z r (n/2 - iz))dz - J kig(Z T (N - iz))dz. 

Now, |Arg(Zr(iV — iz))\ is bounded independently of z when is large enough, thus 
the estimate (3.8) gives the result when 8 > 1. When 5 < 1, we take A^ ~ clog(T) 
with c > 4/£(7 ) and use that log \Z r (n/2 + u — iz)\ < Cz 5 in u £ [0, u }. By what we 
explained above, if Uq is large enough (and we take it independent of T), we also have 
for all z £ R and u > Uq 

| \og{Z r {n/2 + u-iz))\< Ce- u ^ o)/2 

thus 

pN pN 

/ log \Z r (n/2 + u-i)\du+\ log \Z r (n/2 + u- iT)\\du < CT 5 + / e- ue ^ o)/2 du 

Jo Ju 



< CT 



■6 
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and 



J* |Arg(Z r (iV - iz))\dz < ^ e~ m ^ /2 dt < CT' 1 



□ 



Gathering all these results, we obtain (we assume 5 > n/2 as otherwise one has a 
complete formula by (3.7)) 

Theorem 5. Let T be a convex co-compact Schottky group of orientation preserving 
isometries of H n+1 with n + 1 is even and assume that 5, the Hausdorff dimension 
of the limit set, is larger or equal to n/2. Then the Krein function for X = r\EI n+1 
satisfies the Breit-Wigner approximation: for o > fixed and all z G [T — a,T + cr] 
with T large, 

d z i{z) = X (X)P n (z) + - J^% + 0(T 5 ) (3.10) 

where lZ a+ i is the set of resonances in the disc of radius a + 1 centered at T, and P n (z) 
is an odd polynomial of degree n, x(X) the Euler characteristic of X. Moreover the 
following asymptotics holds 

rT 



az)=x(X)P n+1 (z) + R(z), \R(z)\ = 0(z d ), I R(z)dz = 0(T d 
where P n+ i(z) is an even polynomial of degree n + l. 



o 



3.3. The case of surfaces. A convex co-compact surface M = T\M. 2 decomposes 
into a compact hyperbolic surface with totally geodesic boundary (the convex core N) 
and a collection of nf ends which are funnels Fj isometric to 

[0, oo) t x (R/£jZ)e with metric dt 2 + cosh(t) 2 ^ 2 . 

This can be seen as half of the hyperbolic cylinder (7)\H 2 where (7) is the cyclic 
elementary group generated by one single hyperbolic isometry with translation length 
£j. Guillope-Zworski [GuZw] studied the spectral shift function s(z) , related to the 
pair of operators (Am,Af) where := ©"^A^ where Ap. is the Laplacian on 
Fj with Dirichlet condition at the boundary of Fj (given by t = in the coordinate 
above). They define s(z) = ^- log det 5(1/2 + iz) for z £ R, where S(s) is the relative 
scattering operator, and they show that s'(z) is meromorphic in z G C with poles and 
zeros given essentially by resonances of the Laplacians. We do not give details and 
refer the interested reader to [GuZw]. In fact, this function s(z) is given by 

s'(z)=Z'(z)-e F (z) (3.11) 

where £ is the spectral function of (3.4) for M and £f = J2j %Fj{z) where £p. is defined 
the same way as (3.4) on the funnel Fj (and for the spectral measure of A_p.); see the 



1 In [GuZw], it was denoted cr(z) instead of s(x 
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proof of Proposition 4.4 in [BoJuPc]. It is computed by [BoJuPe, (2.10)] that there is 
a constant c such that 

JMtr.(z)+c = Z Fj {l/2 + iz) 

' Z F .{l/2-iz) 

where Zp.(s) := e~ sij ' 4: rifc>o(^ ~~ e -( s +2fc+i)£^2 j g ^ e Selberg zeta function for the 
funnel Fj. This function converges absolutely for Re(s) > and in particular on the 
line 1/2 + z'R. From this convergence, we directly see that £ Fj (z) = ^(-0 an d we a l so 
deduce from Lemma 3.3 (as we did before for Z r ) that £fAz) — f° r £ £ R- Thus 
by combining with Theorem 5 and the Gauss-Bonnet formula Vol(K) = — 2ttx{M) = 
2n{2g — 2 + rif) (where g is the genus), we obtain 

Theorem 6. The scattering phase defined by Guillope-Zworski [GuZw] for hyperbolic 
convex co-compact surfaces M = T\M. 2 satisfies the asymptotics as z — >■ +oo 

s(z) = —Vo\(K)z 2 + R{z), \R(z) \ = 0(z s ), [ R(t)dt = 0(z s ) 
4?r Jo 

where K is the convex core of M , 5 is the Hausdorff dimension of the limit set of V . 
Moreover the Breit-Wigner formula holds: for all z G [T — a, T + a] with T large and 
a fixed 

Air 7r ' \z — or 

In [GuZw], an asymptotics with a remainder O(z) using the method of Melrose [Me88] 
was proved. Our remainder estimate is thus always better than [GuZw] since 5 < 1. 

3.4. Lower bounds for the remainder. For arithmetic hyperbolic surfaces, one can 
use the argument of Jakobson-Naud [JaNa] to get a polynomial Omega lower bound 
on the argument of Selberg zeta function on the line Re(s) = 1/2: 

Proposition 3.6. [Jakobson— Naud] Let F be a convex co-compact subgroup ofisome- 
tries ofM 2 , which is an infinite index subgroup of an arithmetic group Tq derived from 
a quaternion algebra. Assume that the Hausdorff dimension 5 of the limit set satisfies 
5 > 3/4, then for all e > 0, the argument Arg(Zr(l/2 + iz)) of Selberg zeta function 
for z el is not an 0(z 28 ~ 3/2 ~ £ ). 

Proof. Let Zr{\) be the Selberg zeta function of T. In [JaNa, Proposition 3.4], it is 
shown that there exist some constants A > 0, C > such that for all T > large and 
a = 2 log(T) - A, then 

3T QjiiS-2 

\S a , t \ 2 dt>— — (3.12) 
(logT) 2 
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where S att is defined for a > 0, t > as follows: let <p G C£°((— 1, 1)) non-negative 
equal to 1 in [—1/2, 1/2], then set <^ i)Q (x) = e~ ltx ipo(x — a) and 



C((l + \t + z\)~ N ) for all iV G N, we get for t > large 

IS* I < ct 2& - ? ^- s 

where C does not depend on a. This estimate contradicts (3.12) for T large enough, 



Note that such groups exist, see [JaNa, Section 4]. By (3.11) and (3.5), we then 
deduce the 

Corollary 3.7. If V is a convex co-compact infinite index subgroup of an arithmetic 
group T derived from a quaternion algebra, and if the Hausdorff dimension S of the 
limit set ofT satisfies 5 > 3/4, then for all e > the remainder R(z) in the scattering 
phase asymptotics of Theorem 6 is not an 0(z 25 ~ 3 / 2 ~ £ ) . 

In general, without the arithmetic assumption, the argument of Jakobson-Naud [JaNa] 
shows that for all e > 0, if S > 1/2 then the remainder R(z) is not an 0((\ogz)~s e ). 
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